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PREFACE 


The work described in this report was performed by the Mission Analysis 
Division of the Jet Propulsion Laboratory. 
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Abstract 

This report describes the use of a sequential least squares 
filter In the orbit determination for the Mariner Venus-Hercury 
(Mariner 10) spacecraft. The orbit determination strategy outlining 
the use of both the sequential filter and a conventional batch filter 
is given. Highlighted are the mission events from launch to the first 
Mercury encounter with emphasis on the sequential '*ter performance. 
Advantages tc the mission derived from the seque’ .ial filter are 
pointed out. 
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Performance of the Square Root Information Filter 
for Navigation of the Mariner 10 Spacecraft 

The flight of Mariner 10 was launched in November 1973. It 
encountered the planet Venus in February 1974, and the gravity assist 
of that flyby propelled the spacecraft towards Mercury. It encountered 
Mercury in late March 1974, and due to the 2:1 ratio of the spacecraft 
heliocentric period to that of Mercury, it was able to twice again encounter 
Mercury in September 1974 and in March 1975. At each of these four 
planetary encounters a wealth of scientific data was produced including 
thousands of closeup television pictures of the tiny planet Mercury. 

One important ingredient in the success of this mission, the first 
spacecraft to achieve multiple planetary encounters, was highly accurate 
navigation. In performing the navigation, extensive use was made of a batch 
sequential filter and smoother, again, a first for interplanetary spacecraft 
navigation. This report will concentrate on the use of the batch sequential 
filter during the flight, the strategies followed, strengths and weaknesses 
of the filter, comparisons of the filter with a batch filter, and specific 
contributions to the mission from the filter. 

Section I highlights characteristics of the mission and the 
premission analysis that led to Che incorporation Oi the sequential filter 
into Che mission software. 

Section II gives the orbit determination strategy followed in using 
both a classical least squares constant parameters batch filter and the 
new sequential filter which is based on a process noise model; data sets, 
solutions sets, and choices of filter parameters are discussed. 

Section III is a chronology of the mission navigation through the first 
Mercury encounter, highlighting the uses and advantages of the sequential 
filter. 
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Finally, In Section IV Is a discussion of the payoffs derived from 
use of the sequential filter and some thoughts on Its future use. 

The mathematical basis for the design of the filter Is given In 
Reference 1. The specific filter and smoother equations used In the 
Mariner 10 software are given In the Appendix. 

I. Premission Analysis and Planning 

An accuracy analysis study was conducted prior to the mission In 
order to determine the accurate capabilities of the orbit determination 
(O.D.) system, (see Reference 3). This study concentrated on the prime 
mission, l.e., up to the first Mercury encounter. The capability of 
subsequent Mercury encounters was considered, however. 

It was clear from the analysis that the most critical period for 
orbit determination was the period determining the Venus encounter 
conditions. The reason for concern was the multiplication factor of 
approximately 1000 that existed between the Venus and Mercu:ry aiming 
zone errors. This large factor was caused by the large gravity bending 
occurring at Venus encounter. This meant that a 1 kilometer error at 
Venus, if uncorrected, would result in about a 1000 kilometer error at 
Mercury. 

Four trajectory correction maneuvers (TCM's) were planned: TCMl 

at ten days after launch, to remove launch errors; TCM2 at 18 days prior 
to Venus encounter, to control the Venus encounter conditions; TCM3 at four 
days after Venus encounter, to remove Ven. encounter errors; and TCM4 at 
24 days after Venus encounter to remove errors introduced by the expected 
large TCM3 and to control the Mercury delivery. The orbit determination 
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prior to TCM2 was critical in assuring the accuracy of the Venus encounter. 

Table 1 presents lo error ellipses in the Venus encounter plane 
derived from the O.D. during the long arc prior to TCM2. The table 
is composed of the major error sources and the associated 'osition errors. 
Data noise represents the data quality of the doppler and range data*. 

The doppler noise is conservatively considered at 0.015 hz. (lo) or 
equivalently 1 mm/ sec for a 60 sec. count time, and range is considered 
at 50 meters (lo) . Ephemeris error, here, is the effect of the Venus 
position uncertainty relative to the earth at encounter time. Effective 
Station Location Errors (ESLE ) includes uncertainties due to three factors. 

1) The uncertainty of station locations on the earth's crust, 

2) The uncertainty in the position of the earth's crust with 
respect to the earth' s center (earth wobble) , and 

3) The uncertainty due to transmission effects - both in the earth's 
troposphere and ionosphere and charged particles in space (space 
plasma) . 

The troposphere changes the signal path length, and the charged particles 
change the group and phase velocity of the electromagnetic wave. All of 
these disturbances appear to, in effect, move the location of the tracking 
station. 

The dominant error source as shown in Table I is the Non>Gravitational 

Acceleration ( NGA) of the spacecraft. These accelerations are very small, 

_9 2 

typically on the order of 10 m/sec . Accelerations of this size would 
physically perturb the trajectcry by only 3.5 kilomstcts in 30 days. 

^Doppler and range data are defined and discussed in Section II. 
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However, the random time variations of the accelerations make their 
determination a difficult problem at best* The attitude control valves 
(2 for each axis) fire at various rates and leakage can be caused by 
improper valve seating due to dust particles, etc. Additionally, any 
unmodeled spacecraft configuration change can cause a change in solar 
pressure. The error shown in Table I is the result of modeling the NGA's 
as stochastic forces* and using a batch filter (which incorrectly assumes 
all parameters are constant over the entire data arc). 

Table II presents the same information for the Mercury encounter. 

Again note the dominance of the nongravitational accelerations. The liming 
zone or 'science success zone' at Mr .cury was about 500 x 200 km., so the 
predicted performance would meet that objective. Another reason for minimizing 
the Mercury flyby error was to fly close to the 'free return contour'** and 
minimize the fuel needed for a second Mercury encounter. 

As a result of the preflight analysis leading to results such as 
those shown in Table I, the propellent tank size was doubled from a capability 

of 60 m/sec to a capability of 120 m/sec. Still, since about 1/3 m/sec 

would be needed at TCM to remove the error induced by each kilometer of miss 
at Venus, a 3o error at Venus would have necessitated burning all of the 
propellent to correct the error induced and achieve the desired Mercury aim 
point. There would be no chance for a secoiid Mercury encounter. 

* 

The accolerations were modelled as constant over one day, exponentially 
correlated over time with a 3 day correlation time. The a priori uncertainties 
were: 0.5 x 10*^, 1.75 x 10”^, 3.0 x 10" m/sec^ for the three axes. 

"kit 

The free return contour was a locus of points In the Mercury aim plane such 
that an encounter passing through one of these points would have, after 
being perturbed by the planet, a heliocentric period exactly twice that 
of Mercury. 
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The batch sequential filter was designed and Implemented In order 
to better model the NGA^s and reduce the expected errors in the Venus 
delivery. The filter formulation makes use of the square root information 
estimation technology developed at the Jet Propulsion Laboratory and is 
described In the Appendix. 


II. Orbit Determination Strate gy 

The primary data type used for outer space navigation is two-way 
doppler. A radio signal Is transmitted from a station on the earth to the 
spacecraft and retransmitted to the same station (hence, two-way). Doppler 
data consists of a count of the difference between the number of cycles 
received and the number transmitted over a specified count interval (nominally, 
1 minute) • This gives a measure of the average topocentrlc (station centered) 
range rate over the count Interval. Except for periods when the spacecraft 
is near a planet where gravitational acceleration is high, the doppler data can 
be averaged over a longer time (10 to 20 minutes) and still contain all the 
information necessary for orbit determination. This allows the information 
to be compressed into fewer data points. 

A second data type, topocentrlc range, is obtained by measuring the 
time delay for a signal to travel round trip from the station to the space- 
craft and return. Range data is complementary to doppler data in the sense 
that it has different information about the trajectory. Whereas doppler 
data was continuous, at the sample rates mentioned above, throughout the 
mission, range data was more sparse (typically one to three points during 
a station track of 8 - 1C hours, and there were many days with no range data). 



JPL Technical Memorandum 33-757 


5 



Solutions were obtained using three data sets: doppler only, doppler 

and range, and range only, (kinslstency of the solutions based upon the 
different sets Increased confidence In both the data end the solutions. 

Different data arc lengths are sensitive to different error sources. 

ESLE effects and short term spacecraft leaks are much more Important for 
short arcs (2-3 weeks or less). Long arcs are Insensitive to ESLE effects, 
but can be degraded by long term spacecraft nongravltatlonal forces. Hence 
whenever possible, solutions were obtained using data from both the longest 
arc available and a short arc (typically 20 days) containing the most recent 
data. By the longest arc available Is meant data from the most recent large 
crajectory perturbation (e.g., a TCM) to the current time. 

Orbit determination solutions were predicated upon a solution set based 
upon filter models with differing lists of parameters. This was done in 
order to assess the contributions of the various error sources as well as 
the continuing adequacy of the models. A basic consistency among the various 
solutions assured the latter. A solution set consisted of seven solutions 
estimating In turn. 

a) The spacecraft state, 

b) state + station locations, 

c) state + solar pressure, 

d) state + attitude control forces, 

e) state + station location + solar pressure, 

f) same as e) + attitude control, 

g) all parameters. 

Depending on the proximity to a planet, planetary or earth-moon mass and 

planetary and earth-moon epheme:rides were Included In g) . 
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The solution set was run with two filters ; the batch filter and the 
sequential filter. Both filters were mechanized using the numerically stable 
and accurate Householder orthogonal transformations. The mechanization is 
termed the Square Root Information Filter (SRIF) . This mechanization 
is discussed in Reference 1 and more fully in Reference 5. The batch 
filter simply included all the data in a single batch. The sequential filter 
had a process noise model and was able to process sequentially at a fixed small 
batch size (cf Appendix) . 

The sequential filter was run with the three axis attitude control 

accelerations as colored noise stochastic parameters. A batch size of one 

day was used normally and correlation times of from one to five days were 

assumed. (The correlation tlica variation seemed to have little effect on the 

-92 -82 

performance.) A priori sigmas of 10 m/sec and 10 m/sec were used. 

Planetary targeting is usually expressed in terms of the B-plane. 

The B-plane (Reference 3) is a plane passing through the center of a target 
planet perpendicular to the incoming asymptote of a spacecraft. B*T is the 
intersection of the B-plane with the ecliptic and B*R is a vector in the B-plane 
perpendicular to B*T making a right handed system R, S, T where S is the 
incoming asymptote. 


III. The Flight from Earth to Venus to Mercury 

The flight path from Earth to Venus and to the first Mercury encounter is 
shown in Figure 1. Also shown if ':he location in the flight path of the three 
trajectory cC'^rection maneuvers. 

When the spacecraft was within the strong gravity field of the Earth, 


its orbit was determined rather quickly compared to Che time required in 
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interplanetary space. One day after launch when tracking stations around the 
Earth had supplied tracking data, a reasonable 0. D. was obtained. This fit 
was updated over the next few days and used to plan lOMl at launch plus 10 days. 

Using doppler and range data, the one sigma error ellipse In the 
Venus B-plane was about 400 km by 200 km. This was small compared to the 
TCKL execution one sigma error of 1500 km. Tl\e batch filter and sequential 
filter performance was similar in this phase due to the strong near Earth 
data and the short time span. 

The mission phase from TCHl to TCM2 was the most critical for orbit 
determination. The accuracy of the orbit determined the parameters for 
TCM2 which in turn determined the Venus delivery accuracy. Due to the 
geometrv of the swlngby trajectory there was a 1 to 1000 mapping of errors 
from Venus to Mercury. In other words, a one kilometer miss at Venus mapped 
into a one thousand kilometer miss at Mercury. 

A time history of solutions in the Venus B-plane is shown in Figure 2 
for batch filter solutions, and in Figure 3 for sequential filter solutions. 
Solutions for both long and short data arcs are shown for doppler only and 
doppler + range. The long data arc is from TCKI (11/13/73) to ttie date 
shown on the abscissa. The short data arc Is the 20 days preceding the date 
shown on the abscissa. All solutions shown Include the spacecraft state, 
solar pressure parameters, and station locations. This solution set was 
chosen as representative; other solution sets show similar behavior. Note 
that all of the fits with data up to 12/24 showed good agreement. The batch 
filter fits with data up to 1 ^ 2 , however, show a wide scatter in B*R. One 
week later with data up to 1/10 the wide scatter is apparent in both B*R and 
a*T with the short arc doppler solution 800 km high in B*R. The sequential 
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filter, on the other hand, with the exception of a short arr solution high 
in B*R, remained remarkably flat throjghou. ;he entire time period. 

Some insight into the scatter of solutions shown in Figure 2 and 
the apparent stability of the sequential filter solutions in Figure 3 
can be gained from examining Figure 4. This figure is a plot of the space- 
craft nongravitational accelerations as represented by the magnitud ; of 
the smoothed solution vector for the attitude control acce^ '^ations. Of 
course, this is simply the acceleration that the filter sees, i.e. , that 
along the Earth-spacecraft line of sight. The filter parameters were: 

one day batch size, 5 day correlation times, and a priori sigmar of 
-9 2 

. X 10 m/sec along each axis. 

Significant peaks in Figure 4 are numbered. Peak Nos- 1, 2 and 3 
occur on days when roll calibrr.tions were performed. The spacecraft was 
allowed to roll about the spacecraft - sun axis at a rate of 40 rain/rev 
for 6 revs. This meant that for 4 hours the solar pressure component 
perpendicular to the sunline was averaged out, so the average acceleration 
due to this component was less for the day. When it was first noticed 
that the sequential filter was sensitive to this minute force, the Navigation 
Team members were surprised and gratified that the filter and the model were 
performing so well. 

Peak No. 4 occurred on January 4. This unusually high acceleration down 
by the filter was confirmed by a ramp in the doppler lesidualF beginning on 
January 4 and lasting for about a day. The total doppler shift was about 
-0.01 hz. corresponding to a velocity change of 0.6 mm/sec. Subsequent analysis 
of telemetered limit cycle data indicated a pitch axis torque lasting for 
about 25 hours on January 4 and 5 . It was calculated that a valve leak 
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large enough to cause the observed torque and the observed velocity change 

would have been less than 0.9 grams of gas e 3 q>ended. This is coo small tc < 

be detected by the gas comsumption telentetry. 

A 

Here is an example of a spacecraft nongravitational force (a gas 
leak) too small to be seen by telemetry but large enough to effect the 
orbH determination (the scatter of solutions shown in Figure 2). The 
sequential filter modelled the disturbance beautifully and provided a 
stable Venus B-plane solution to use in the determination of the parameters 
for TCM2. 

Peaks 5 and 6 have no easy explanation. Peak 5 is probably a filter 
induced damping to fit of peak 4 and peak 6 is probably another sp ^ecraft 
leak. 

After TCM2, 0.1). solutions were made every two days in order to get a rapid 
determination of the trajectory. Had TCM2 not been successful there was 
still time to do another roanerver prior to the Venus encounter. Here the 
batch and sequential filter gave similar results over the one week period. 

Figure 5a shows the solutions in the Venus B-plane for 3, 5 and 7 day arcs. 

On January 28, one week before Venus encounter, the spacecraft 
developed a serious problem. The roll gyros began oscillating at a rapid 
rate, causing the roll position jets to fire at their maximum rates. The 
immediate effect was a large consumption of cold nitrogen used as attitude ^ ' 

control gas (0.54 kg. in 79 minutes) before the problem was diagnosed and the , 

f. 

gyros turned off. The thrust imbalance in the attitude control iets imparted I 

a force on the spacecraft. The doppler pseudoresidualo'^ indicated a | 

1 : 

doppler shift of 0.13 hz. or a radial velocity change of 8.5 mm/sec in 79 minutes, | ' 

} ■ 

— ^ ' 

* i 

observed-computed values 
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hence the radial acceleration was 1.79 x 10 m/sec . This was a large 
eiough acceleration to have two effects: 1) the O.D. data arcs had to be 

restarted, and 2) the trajectory was changed thus altering the Venus 
B-plane position. It was suspected that at the tine that the unseen 
component of the acceleration vector could have been large enough to move 
the trajectory up to 100 km in the Venus B-plane. That size error at 
Venus would have seriously degraded the science return at Hercury. 

After the gas acceleration on January 28, short arc fits were again 
made every two days. Figure 5b shows these solutions for 3, 5 and 7 days. 
Also shown is the actual trajectory point determined after the Venus flyby.* 
By con^aring these solutions with those cf the previous week, it is seen 
that the January 28 acceleration moved the trajectory only about IS km in the 
B-plane. Tills number is consistent with that obtained postflight by solving 
for the acceleration and analyzing the spacecraft telemetry giving gas usage 
and torques. 

Also shown in Figure 5 is the desired aim point and the actual aim 
point. These are different because the spacecraft had discrete turn and 
burn time capabilities. The actual aim point was the closest realizable 
aim point to the desired point. The proximity of the flyby point to the 
desired aim point (17 km) is remarkable. It is much smaller than the 
preflight statistics shown in Table 1 indicated. This was fortuitous in 
that, because of the small miss, a large TCM shortly after Venus encounter 
was not necessary. Further the direction of the miss was such that the 

. to the Venus gravity field experienced approaching and leaving 

planet, the statv of the spacecraft was determined very accurately 
' within less than i. km in the Venus B-plane) . 
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desired Mercury aim point could be achieved with a * sunline** TCM or 
March 16. 

On February 14 the gyros were turned on in order to test the conditions 

under which the January 28 failure had occurred. The oscillations again 

"“6 2 

appeared and another spacecraft acceleration resulted (1.75 x 10 m/sec 

for 20 min.). It was then decided not to attempt a full maneuver but to lise the 

sunline TCM on March 16. 

Figure 6 is a continuation of Figure 4 showing the radial nongravi- 
rational acceleration from TCMl through the Mercury encounter. It was 
constructed by joining 5 different O.D. fits (data arcs starting at each 
TCM and at the large anomalies of January 28 and February 14) . Ti additional 
peaks are noted. Peaks marked 7 and 8 on days 107 and 123 correspond to 
* flyback and sv^eeps* where the reference star was lost and attitude control 
gas was expended in the reacquisition search. It is seen that the level 
of the spacecraft nongravitational accelerations is generally higher after 
the January 4 peak. 

An example of the effect of this high level on the batch filter can be 
seen in Figure 7 which shows various O.D. solutions in the Mercury B-plane. 

The solutions labelled 1 contain data only near Venus encounter (February 4) . 

The solutions labelled 2 and 3 contain data starting after the February 14 
leak (day 103) • The batch filter solutions are badly scattered whereas 


The engine nozzle pointed toward the sun in the normal spacecraft cruise 
attitude, so if the engine is fired with no turns, the thrust vector 
imparted to the spacecraft is along the sun-spacecraft line, hence the 
name 'sunline* • Rather / three degrees of freedom, the sunline maneuver 
only has two: the magnitude of the thrust and the time the maneuver takes place. 

The advantage of a sunline maneuver was that the gyros needed only be turned 
on for the few seconds during which the engine fired, rather than the hour or 
so needed for a full maneuver with turns. 
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the sequential filter solutions cluster relatively wlghtly. Th.i scatter 
of the batch filter solutions can be attributed to ^ts Inability to model 
the high level of nongravltational accelerations. The data spans for 
solutions 2 and 3 are indicated in Figure 6. It Is seen that solution 2 
contains peak 7 and solution 3 contains both peaks 7 and 8. As with TCM2, 
the best orbit used to determine the TCM3 parameters was given by the 
sequential filter. 

After XCM3 the primary purpose of orbit determination was to accurately 
determine the trajectory so that the science Instruments could .. pointed 
accurately. The Mercury flyby was about 160 km from the aim point, trell 
within the science success zone. 

The spacecraft was now In an orbit about the sun %d.th a period twice 
that of Mercury. This enabled two more Mercury encounters, one in 
September 197A and one in March 1975 before the spacecraft attitude control 
gas was depleted. The sequential filter was used as well as the batch 
filter throughout the remainder of the mission. The performance differential 
of the two filters, however, was not as dramatic as that in the early mission 
described above. Some navigation results of the extended mission are given 
in References 4 and 6. 
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IV. Conclusions 

The performance of the square root mechanized batch-sequential filter 
in the navigation of Mariner 10 was impressive. It demonstrated that the 
SRIF formulation is a viable one for use in mission operations and also that 
a first order colored noise model adequately described the NGA effects. 
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The sequential filter was heavily used and relied upon during the 
nd-sslon. It Is difficult to say what would have transpired had It not been 
avs .lable* or a different filter mechanization used bur It Is quite conceivable 
tha a less accurate O.D. during the critical pre-Venus phase would have 
caused a large change In the remainder of the mission. Thou^ the direction 
of t. he error at Venus was fortuitous in that a sunline maneuver was possible 
to obtain the desired Mercury delivery, a large error at Venus %rould have 
severely degraded the Mercury delivery. This, in turn, would hav«. not only 
degraded the science return on the first Mercury encounter, but could have 
precluded any further Mercury encounters. 

Based on .he success Indicated here, the square root sequential filter 
\;itl process noise is now a tested, accepted tool for interplanetary orbit 
determination. Essentially the same . .Iter with some added capabilities 
(e.g., variable batch sizes, ability to model white noise as well as 
colored noise, srootned covariances, and smoothed sensitivities will be used 
In the upco’<*ng Viking mission to Mars. This state-of-the-art estimation 
mechanization Is described in Reference 1. 
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Table I 

Venus Preencounter Error Budget 


I Sigma Error Ellipse 

Error Source in the B-plane (km) 

SMAA SMIA 


1. 

Data Noise 

18 

X 

2 

2. 

Ephemeris 

9 

X 

5 

■» 

^ • 

ESLE 

38 

X 

5 

4. 

NGA 

93 

X 

22 

5. 

Total 

99 

X 

36 


(SMAA - semimajor axis, SMIA - semiroinor axis) 


Table II 

Mercury Preencounter Error Budget 

1 Sigma Error Ellipse 
Error Source in the B-p)ane (km) 


1. 

Data Noise 

SMAA 

10 

X 

SMIA 

7 

2. 

Ephemeris 

65 

X 

42 

3. 

ESLE 

20 

X 

8 

4. 

NGA 

159 

X 

120 

5. 

Total 

172 

X 

129 
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Fig. 1. Projection of the Earth to Mer- 
cury trajectory into the ecliptic plane 
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Fig. 2. Time history of solutions in the 
Venus B- plane: batch filter g 


Fig. 3. Time history of solutions in the 
Venus B-plane: sequential filter ^ 
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Appendix 


ential Filter Formu 
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In practice M was the diagonal matrix, M - e I, where At is the batch 
. ^ . -4 T -to •lapntitv matrix. Hence the 


size, T the correlation time, and I is the identity matrix. Hence the p 
stochastic parameters! are uncorrelated from each other but each is correlated 

in time. 
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Equations (A. 5) and (A. 6) form a backwards recursion for x* and p* 
that. If started with the final filtered estimates (which are, by definition, 
smoothed estimates), yield smoothed estimates for each batch. 
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